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Abstract 

We present a comparison of our model treating fission dynamics in conjunction with 
light-particle (n,p,a) evaporation with the available experimental data for the nuclei 
126 Ba, 188 Pt and three isotopes of the element Z=110. The dynamics of the symmetric 
fission process is described through the solution of a classical Langevin equation for a 
single collective variable characterizing the nuclear deformation along the fission path. A 
microscopic approach is used to evaluate the emission rates for pre- fission light particles. 
Entrance-channel effects are taken into account by generating an initial spin distribution of 
the compound nucleus formed by the fusion of two deformed nuclei with different relative 
orientations. 
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1 Introduction 



The dynamical time evolution of the fission process from an initially formed compound nucleus 
with a more or less compact shape to the saddle and scission configurations and the simultaneous 
emission during this deformation process of light particles constitutes a very complex problem. 
It is the aim of our study to describe this process taking into account the excitation of both 
thermal and rotational nature, the initial configuration and to calculate the evaporation of light 
particles from an excited, rotating, deformed nucleus on its way from the initial compact shape 
to the scission point. 

In the absence of a complete microscopic ab initio theory of such a dynamical process, many 
different theoretical approaches aiming at its description have been worked out [§J-0. They 
generally rely on a classical description of the evolution of the collective coordinates which are 
introduced as abundant variables. In these descriptions, collective parameters appear (collective 
mass, friction, and diffusion coefficients) which depend on the collective coordinates. It seems 
clear that the quality of the theoretical description will depend on the more or less pertinent 
choice of the collective coordinates and the degree of realism of the underlying theory used to 
determine the collective parameter functions. The importance of pure quantum effects such 
as pairing correlations and the existence of shell effects which are present at low excitation 
energies has been investigated recently ]T0|]-|p|. 

The multiplicity and the characteristics (energies, angular distributions, correlation func- 
tions) of prefission light particles and gamma emission may provide some useful information on 
the time evolution of the nucleus as it evolves towards the saddle and scission point. Observables 
like prefission light-particle multiplicities may, indeed, work as a clock. Angular distributions 
and particle correlations may give informations on the surface and deformation of the fissioning 
nucleus. 

We recently developed a realistic time dependent model which described the time evolution 
of an excited nucleus as it is e.g. generated in a heavy-ion collision and decaying through sym- 
metric fission with pre- and post-fission light particle emission JTSJ. We applied our description 
to the case of 160 Yb and made a detailed analysis of different physical quantities characterizing 
the decaying system. We found a rather good agreement of our theoretical predictions with the 
experimentally observed light-particle multiplicities. 

More recently, several new sets of experimental data have become available 16fl-[21| with 
a more complete analysis of the emitted light particles. It is the aim of the present work to 
systematically test our approach by a comparison with these experimental data. The fact that 
these data are obtained in quite different regions of the periodic table makes a comparison 
between calculated and experimental results all the more challenging. 

A particular motivation for this work is the explicit description of the entrance channel 
through which the initial spin distribution of the system is generated. This is achieved by 
considering all possible relative orientations of the two colliding deformed nuclei in the entrance 
channel. Their different orientations and the different possible impact parameters of the reaction 
yield the initial spin distribution da/dL versus L which determines the relative weights of the 
initial angular momentum of the compound nucleus with which the Langevin trajectories are 
started. 

Still another motivation for the present work is the presentation of the link between the 
well known Weisskopf formula |E| and a new and more microscopically founded description of 
the emission width for the light particles in terms of phase-space densities of the particle to be 



emitted |23| 
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The paper is organized as follows: The dynamical model describing the time evolution of 
a deformed, hot, and rotating nucleus from its initial to a final configuration is presented in 
section 2. In section 3 we specify how to determine the emission rates for light particles from a 
hot, deformed and rotating nucleus moving along the previously described trajectory. In section 
4 we show how we obtain the initial spin distribution which has been not considered in our 
previous study ||15 |. Our results concerning the nuclear systems 126 Ba, 188 Pt and 266 > 272 > 278 no 
are presented in section 5. The paper is closed in section 6 with a summary of the most relevant 
results and an outlook on planned extensions of the model. 



2 Collective dynamics of an excited system 



We consider an ensemble of deformed nuclei with finite excitation energies and rotational an- 
gular momenta as given by the initial conditions determined from the entrance channel. The 
subsequent time evolution of the nucleus is governed in our present description by a single 
collective coordinate q = Ri 2 /Ro where Ri 2 is the distance between the two centers of mass 
of the left-right symmetric deformed nucleus and Rq is the radius of the corresponding spher- 
ical nucleus having the same volume. This collective variable is defined in the framework of a 
Trentalange-Koonin-Sierk (TKS) ]24| parameterization of the surface of the nucleus. The TKS 
deformation parameters are related to q by means of a minimization procedure of the collective 



potential energy defined below [25 



Denoting the conjugate momentum by p(t) we use the following classical equations of motion 
to describe the time evolution of the fissioning nucleus 
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dM{q) dV{q) 7 (g) 



dq 



dq M(q) 



P + F L (t) 



(1) 
(2) 



Here M(q) is the collective mass determined in the incompressible fluid approximation [g(J an d 
7(g) the friction coefficient calculated in the wall-and-window friction model ||27| . The collective 
potential V(q) is obtained as the difference of the Helmholtz free energy at deformation q minus 
the one for the ground-state deformation. 

The free energy could in principle be obtained from a microscopic mean-field calculation at 
finite temperature. To perform such a constraint Hartree-Fock calculation using a reasonable 
effective nucleon-nucleon interaction of the Skyrme or Gogny type at every point in the multi- 
dimensional deformation space is, of course, completely out of question due to the tremendous 
computer time such an analysis would involve. Even to perform the same kind of calculation 
on the level of a selfconsistent semiclassical approximation like the Extended Thomas-Fermi 
(ETF) method [^] at finite temperature pSfl , which would describe the average nuclear struc- 
ture without shell oscillations, would be far too time consuming. This is why we have rather 
used a still simpler semiclassical approach and the liquid drop model with the parametrization 
of Myers and Swiatecki |30| as explained in Appendix A. We have tested the barrier heights 
obtained in this way versus the ones obtained from selfconsistent semiclassical calculations (at 
zero temperature) and obtained agreement within a few MeV as described in Ref. ||25||. These 



calculations describe, of course, deformation properties of non excited nuclei. In order to take 
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nuclear excitation into account we have used the temperature dependence of the LDM parame- 
ters associated with the Skyrme SkM* interaction |31[|, determined for that interaction through 



selfconsistent semiclassical calculations at finite temperature in Ref. |25[ . 

To use the temperature dependence of the LDM parameters associated with the Skyrme 
SkM* force together with the Myers-Swiatecki parametrization of the LDM is, of course, in- 
consistent. Since approaches give very similar results for the semiclassical energy at zero tem- 
perature this approximation seems, however, very reasonable. The semiclassical approach used 
here is of course only an approximation at low temperatures since shell effects are absent from 
our description. It is however well known that nuclear shell effects are washed out with in- 
creasing nuclear temperature and have essentially disappeared beyond T~ 2.5-3 MeV. At these 
temperatures the semiclasssical results are becoming exact. For the cold systems the fission 



barriers obtained with this potential are higher that those evaluated by Sierk [33], but already 



at excitation energies E* ~ 90 MeV both barriers become comparable and for E*>90 MeV the 
barriers evaluated in our model are even smaller than the Yukawa folded ones of Sierk. 

The friction term and the Langevin force F^t) in Eq. (0) generate the irreversible pro- 
duction of heat energy and the energy fluctuations respectively which both originate from the 
coupling of the collective dynamics to the intrinsic degrees of freedom. In practice one de- 



fines = y D(q)fi(t), where D(q) is the diffusion coefficient. We take the simplifying 

point of view that it is related to the friction coefficient 7(g) through the Einstein relation 
D(q) = 7(g) T, where T is the temperature of the system. The quantity /l(t) can be written 
in the form /l(t) = y/rrj, where r is a time step length corresponding to a time interval [t, t + r] 
and 77 is a gaussian distributed random number with zero average (77) = and variance (ry 2 ) = 2 
where brackets represent ensemble averages. 

The Einstein relation is in principle valid at high temperature, in a regime where the process 
can be described in a classical framework. Quantum effects (pairing, shell corrections, collective 
shape vibrations) may be present at low temperatures and modify the relation between D and 
7 [lll]-||14||. We shall come back to this point later on. 

In principle, one would have to treat each ensemble of nuclei with a given initial angular 
momentum L and a given initial excitation energy microcanonically. We assume that we may 
replace this microcanonical ensemble by a grand canonical one in which, instead of the mean 
excitation energy, the temperature T has a well-defined value at a given time. This simplifying 
assumption is innocuous as far as the mean values of observables are concerned, but it may 
falsify the fluctuations to a certain extent. We assume that the time scale which governs the 
fission dynamics is much larger than the internal equilibration time, otherwise the definition 
of a nuclear temperature would make no sense. Under these conditions the system can be 
considered as being continuously close to equilibrium. We suppose that the nuclear excitation 
energy E* is related to T through the usual Fermi-gas expressions E* = a(q)T 2 , where a(q) is 
the level density parameter at a given deformation q. 

For practical calculations Eqs. (0) and (0) are rewritten in a discretized form and numeri- 
cally integrated over small time steps r |Tj| . In order to perform such an integration one needs 
to start from some fixed initial conditions which define the beginning of the process. We take 
into account the possible emission of light particles in every time step and we make sure that 
the average total energy of the system will be conserved. This determines the nuclear excitation 
energy and hence the temperature at each instant of time. 
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3 Particle emission from a hot, deformed, and rotating 
nucleus 

Particle emission before scission is governed by transition rates T^(E*,L) which determine 
the number of particles of type v (we take into account neutrons, protons, and a particles) 



+ \Ae 



and with an 



emitted per unit time with an energy e a in an interval 
angular momentum £p from a nucleus with average excitation energy E* and total angular 
momentum L. In Ref. [I5| we used the well known Weisskopf formula \T1\ for the partial 
width T"P(E*,L) in terms of densities of states of the emitting and residual nucleus and of 
the transmission coefficient w v (e,£p) for emitted a particle v with given energy e and angular 
momentum £p. The determination of w u takes into account the deformation and rotation of 



the emitting nucleus (see Appendix B of Ref. |L5| for details). 



In Ref. another more microscopic determination of these transition rates was proposed. 
In this framework the transition rates are given as 



d 2 



-Ae A£ . 



(3) 



de a d£/3 

where e a and £p characterize an emission energy and angular momentum lying in the intervals 

and 



2 £ * 



'-As 
2 



'/3 



-A£, £ + -A£ 
2 ' p 2 



respectively. 

The number n v of particles of type v which are emitted per time unit through the surface 
£ of the fissioning nucleus is given by p3f 



(4) 



n v = da d 3 p f I/ (r ,p')v' ± (r' )w l/ (v , ± (r t )) 



where p', v' are the momentum and velocity in the body-fixed frame. The quantity v' ± is 
the velocity component perpendicular to the emission surface at the surface point fg. The 
p' = mv' + mujxr' is the momentum of the particle of mass m in the laboratory reference 
frame and u the angular velocity of the nucleus in this frame. Here and henceforward, primed 
quantities refer to the body-fixed frame. 

The classical distribution in phase space reads 



A 



(f',p') 



cxp 



T \2m T U 



uj£' 



(5) 



The 9 function is 1 if r ' lies inside the nuclear volume or on its surface £ and zero otherwise. The 
quantity fi u is the chemical potential and £' the body-fixed angular momentum in the direction 
perpendicular to the axis of rotational symmetry of the deformed nucleus. The potential U is 
taken as 

U{r') = -V + V Cb (f) , (6) 

where Vq > is chosen as a constant mean field potential and Vcb is the Coulomb potential 
experienced at f ' by protons. 
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The quantity w u (v' ± (r' )) is the classical transmission coefficient for the emission of a particle 
of type v. The transmission factor w v was chosen to be the one of an inverted harmonic oscillator 

The explicit relation of (3) with the Weisskopf formulation is discussed in Appendix B for a 
spherical emitting nucleus with angular momentum L = 0. The rates given by (3) can be worked 
out numerically if the distribution function in the phase space is known, what is the case for the 
neutrons and protons but not for the a-particles. We are working on a model which describes the 
distribution of a-particles by folding of the distributions of four (2n and 2p) correlated particles. 
The present approach is like the one of Ref. [lTSJ and it allows us, in principle, to determine 
the particle emission in a given direction of space, hence the determination of observables 
like angular distributions and particle-particle correlation functions. Such observables may be 
worked out in the future. We did not attempt to do it here because of the lack of corresponding 
experimental data. 

The transition rates Y^f are used in a simulation algorithm by means of which we determine 
at each time step [t, t+r] along each classical trajectory whether a particle of given type with an 
energy and angular momentum in given intervals is emitted or not from the compound nucleus. 



Since the algorithm is already described in detail in Ref. [15| we do not repeat it here. 



3.1 Initial conditions and energy balance 

In order to integrate the classical equations (1) and (2) we need to fix the initial conditions from 
which the compound system starts and evolves either through a fission channel to its saddle 
and scission point or stays as a compound system which only emits light particles, i.e. ends up 
as an evaporation residue. 

All the experimental systems which are considered in the next section are generated by 
means of heavy ion collisions at some bombarding energy. The nuclei which are involved in the 
reaction process can be deformed. The initial conditions corresponding to the origin of time are 
fixed by go and Po, the initial value of the collective variable and its conjugate momentum and 
the spin distribution of the system which fixes the relative weight of the angular momentum of 
the initial compound systems. 

The coordinate go is fixed at the value of g where the collective potential V(q) is minimal 
and its conjugate momentum is drawn from a normalized gaussian distribution 

P(po) = (2nMT y/ 2 exp(-p 2 /2MT ) , (7) 

where M = M(q ) is the collective mass. The initial temperature T is obtained through Eq = 
a(go)TQ , where Eq is the initial excitation energy which can be obtained from the knowledge 
of the total energy as explained below. 

In the reaction process, the compound nucleus can be formed with different values of the an- 
gular momentum. If both nuclei are spherical it is easy to construct the initial spin distribution 
under the assumption that the reaction cross section is given by the geometrical expression 

where k is the wavenumber of the relative motion. 

The situation is more complicated if one or both initial nuclei are deformed. In principle 
one should then consider all possible relative orientations of the nuclei and follow their relative 
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trajectories from an infinite distance up to fusion. This implies the knowledge of the nuclear and 
Coulomb potential for given relative distance between the centers of mass and given orientation. 



This raises intricate problems which have been considered by several authors |34], |36[ but 
not yet solved in their full generality. 

In a first step we want to avoid such cumbersome calculations and restrict ourselves to a 



simplified procedure which has been already proposed in former work [[36], |37|. In practice we 
consider the scattering of spherical ions of fixed energy E and introduce the gaussian energy 
distribution with the width 

AE = 5E + AxE + A 2 E , (9) 

where SE is the experimental beam width, A\E describes the slowing down of the projectile in 
the target, and the quantity A 2 E is the difference of the Coulomb barriers between two extreme 
geometries of touching spheroidally deformed nuclei, i.e. with aligned (tip on tip) or parallel 
(side on side) z-axes. 

One can now work out the classical trajectories describing the relative motion of the equiva- 
lent spherical ions for every initial angular momentum L in [0, L max \, where L max is the maximal 
angular momentum for the colliding ions and for a sufficient number of energy values within 
the distribution defined above. 

Repeating the trajectory calculations leads to the spin distribution 

where Lj is the considered angular momentum, Nf is the number of trajectories which lead to 
fusion and iVj is the total number of trajectories. The quantity k is the wave-number of the 
relative motion of the incident nuclei. The present procedure, even though it is not rigorous, 
leads to spin distributions which are rather realistic |36], j37fl . 

Finally, in order to describe the evolution of the excited compound system one needs to 
follow the evolution of the average excitation energy along the trajectory. This is achieved by 
the requirement of energy conservation. 

At the initial point of a trajectory the total available energy can be written as 

Etot = -Ecoll + E ro t + E 

= i^H + v (<io) + ^^ + e* . (ii) 

2M(q ) 2J(g ) 

Here, J(qo) is the moment of inertia of the compound system taken as a rigid deformed ro- 
tator and Eq is the initial intrinsic excitation energy of the compound nucleus. The intrinsic 
excitation energy E*(t) at any given later time t can be determined from the energy balance 

Etot = -E'coii + E Iot + E£ + B u + e a + E Icco n , (12) 

where B u is the binding energy of the emitted particle, different from zero only for a-particles, 
e a is the kinetic energy of the emitted particle and -Brecon is the recoil energy of the nucleus 
after the emission of a particle, which can be neglected in practical calculations. 

For each choice of the initial conditions one generates a separate trajectory. Emitted parti- 
cles are counted with their energy and angular momentum. If the system overcomes the fission 
barrier the trajectory (event) contributes to the final fission cross section 



daa ss ^ / da F \ N 



rfiss 



afiss = ^^Z7 = MrfZ" J "W"' (13) 
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where N^ ss is the number of trajectories which lead to fission, and Nf - the number of fused 
trajectories at a given angular momentum Lj. The sum runs over all angular momentum bins 
and dap/dL is the fusion cross section given by Eq. (10). 

The numbers of prefission particles obtained for each angular momentum of the compound 
nucleus (M v ) are weighed with the differential fission cross section in order to obtain the 
measured number of particles emitted in coincidence with fission: 

<Mu > = — . (14) 

Ofiss 



4 Application of the model to various experimentally 
studied systems 

In the present section we aim to confront experimental results concerning the multiplicities of 
emitted particles in coincidence with symmetric fission events with the model described in the 
preceding section. All measured and calculated data are presented in Table 1 for three systems: 
126 Ba, 188 Pt, and 266 > 272 > 278 no obtained by different entrance channels and at different energies. 
These data result from experiments performed at the SARA (Grenoble) and the VIVITRON 
(Strasbourg) [[U|-[|(J using the DEMON neutron detector |38], j39fl . Unfortunately, the charged 
particles (p, a) were not measured, so we give in Table 1 the theoretical predictions only. We 
present a more detailed discussion of the results for each system in the next subsections. 



4.1 The Ba compound system 

The compound nucleus 126 Ba has been experimentally produced through two different entrance 
channels: 

(I) 28 Si + 98 Mo -> 126 Ba at E lab =20A, 187, and 166 MeV 
and 

(II) 19 F + 107 Ag ^ 126 Ba at E lab =U8 and 128 MeV. 

Fig. 1 shows the fusion cross-section and the fission yields as a function of the initial angular 
momentum. The distributions of the initial angular momentum were calculated with the help 
of the corresponding Langevin equation || || for two different entrance channels (I, II) and 
various bombarding energies. As it can be seen, the fission yields are rather small and located in 
the tails of the distributions. This observation which also holds for many other systems implies 
that it is of great importance to calculate carefully the dependence of the fusion cross-section 
on the angular momentum if one wants to describe the competition between the decay of the 
compound nucleus by fission and by light particle emission in a correct way. Concerning Fig. 1 
we also note that the dominant part of the initial excitation energy resides in the collective 
rotational motion. One observes in Fig. 1 a general increase of the fission cross section with 
the excitation energy of the compound nucleus 126 Ba with the only exception of the highest 
excitation energy. In that case the emission of a-particles becomes important and competes 
with the fission process. This causes a substantial loss of internal energy due to particle emission 
which in turn leads to an increase of the fission barrier and therefore to a smaller fission yield. 

In Fig. 2, the potential energy Va ss of the fissioning nucleus 126 Ba for different angular 
momenta is plotted as a function of the relative distance between fission fragments. One can 
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see in Fig. 2 that the fission barrier of the fused nucleus at temperature T—1.6 MeV and 
angular L=0 is quite high in our model. It is equal to 43 MeV. It vanishes at very large angular 
momentum (L m 80h). The high temperature can also reduce significantly the fission barrier 
as one can learn from Fig. 3, where the free energy of 126 Ba is plotted for a few temperatures 
between and 4.8 MeV. This is why one can observe a significant fission rate only at large 
angular momenta and high excitation energies. All curves representing the fission barriers in 
Figs. 2 and 3 end at the scission point region. It is seen from Figs. 2 and 3 that the saddle 
and scission points are close to each other for 126 Ba. 

Part of the excitation energy of 126 Ba is contained in the rotational mode. In Fig. 4, the 
temperature of 126 Ba corresponding to different excitation energies is plotted as a function of 
angular momentum. One can see that, at the lowest excitation energy E* =84 MeV and the 
highest angular momentum, the temperature is 0.5 MeV only. This means that our theoretical 
predictions based on the statistical equilibrium could be too rough in this case. 

The experimental neutron multiplicities obtained for both the fission and the fusion channel 
are shown in Table 1 along with the calculated neutron, proton, and a-particle multiplicities. 

In Figs. 5 and 6, we show the calculated average numbers of neutrons, protons, and a- 
particles emitted in coincidence with fission as a function of the excitation energy E* of the 
initial compound nucleus. The initial compound nucleus 126 Ba is produced by the fusion reac- 
tion (I) in Fig. 5, and in Fig. 6 the same initial compound nucleus originates from the fusion 
reaction (II). Only the n-emission on the way to fission has been measured so far. The measured 
numbers of prefission neutrons are indicated by points with error bars. At two excitation ener- 
gies we have measurements for the two different entrance channels. As one can read from Table 
1 at the excitation energy E* = 101.4 MeV the measured and calculated numbers of emitted 
prefission neutrons (in coincidence with fission) are M n = 1.32 and 1.83, resp., for the entrance 
channel (I) (Fig. 5) and M n = 1.31 and 1.80, resp., for the entrance channel II (Fig. 6). At 
the higher excitation energy E* = 118.5 MeV the measured and calculated numbers of emitted 
neutrons are M n = 2.01 and 1.71 for the entrance channel (I) (Fig. 5) and M n = 1.85 and 1.99 
for the entrance channel (II) (Fig. 6). 

The experimentally observed number of emitted prefission neutrons at given excitation 
energy E* of the initial compound nucleus is thus larger for the entrance channel (I) than for 
the entrance channel (II) and the difference of the number of emitted neutrons is seen to be 
larger for the larger excitation energy. 

The interpretation of this observation is that for 126 Ba only the compound nuclei formed 
with high angular momentum may undergo fission and thus give rise to prefission neutrons 
because the fission barrier for low angular momentum is high (> 15 MeV). For T—1.6 MeV the 
LD fission barriers vanish for about L = 70 h as can be seen from Fig. 2. Thus the main part 
of the prefission neutrons is obtained for angular momenta around L = 80 h as indicated in 
Fig. 1. At given excitation energy E* the number of compound nuclei formed with such high 
angular momentum is larger for the entrance channel (I) with 28 Si as a projectile than for the 
entrance channel (II) with 19 F as a projectile (see Fig. 1). 

The theoretically estimated prefission proton multiplicities (M p ) are very small as one can 
see in Figs. 5 and 6. The multiplicities of a-particles are not negligible at the highest (E* =131.7 
MeV) and the lowest (E*=8A.l MeV) excitation energies (see Figs. 5 and 6). The increase of M a 
at high excitation energy is mostly due to the high temperature effects while the unexpectedly 
large number of emitted a-particles at E*=8A.l MeV is mostly due to the large deformation of 
the emitter and the huge centrifugal potential (L > 80k) reducing the already small Coulomb 
barrier for the a-particles at the tips even further. 
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As for the comparison of measured and calculated numbers of emitted prefission neutrons, 
the general agreement is not unsatisfactory. Nevertheless, one notices that the calculated 
average emission numbers are larger than the observed ones for the higher excitation energy 
E* = 118,5 MeV. If, instead of the temperature-independent friction, which underlies these 
calculations, we would use a friction parameter increasing with temperature, the agreement 
might be better. 



4.2 The 188 Pt compound system 

As seen in Table 1, the compound nucleus 188 Pt was produced at two excitation energies 
E*=100 MeV and 66.5 MeV using two different reactions 

• 34 S + 154 Sm ^ 188 Pt at E lab = 203 and 160 MeV , 

• 16 + 172 Yb ^ 188 Pt at E lab = 138 MeV. 

The theoretical estimates of the initial spin distributions (solid lines) as well as the fission rates 
(bars) for £'*=100 MeV and two entrance channels as well as for £'*=66.5 MeV are plotted 
in Fig. 7 as functions of angular momentum of the compound nucleus. It is seen that each 
reaction leads to different spin distribution. These differences are responsible for the entrance 
channel effects in the prefission neutron multiplicities emitted by 188 Pt at £'*=100 MeV (see 
Table 1). 

Similarly, as in the case of 126 Ba, we present in Fig. 8 the fission barriers of 188 Pt obtained for 
different angular momenta (L) at fixed temperature T=1.6 MeV. It is seen that for T = 1.6 MeV 
and L = 70 h the fission barrier becomes negligible. Also one can see in Fig. 9 that with 
increasing temperature the fission barrier decreases. At T = 4.8 MeV and L = 30 h, the fission 
barrier of 188 Pt is about four times lower than for T = and L = 30 h. Contrary to the case of 
126 Ba, where the saddle point is very close to the scission point, the path from the saddle point 
(q = Ru/Rq ~ 1.6) to the scission point (q ~ 2.2) of 188 Pt is much longer and will take more 
time. This implies that the case of 188 Pt is better suited to study the influence of dynamical 
effects on the prefission light-particle multiplicities. 

The prescission neutron multiplicities theoretically predicted for 188 Pt are too small (by 
~1 unit on the average) in comparison with the experimental data [TfJ. Also the effect of the 



entrance channel is not fully reproduced. Unfortunately the protons and a-particle multiplicities 
were not measured, so we cannot test the predictive power of our model in this case. 

4.3 The Z = 110, N = 156, 162, and 168 compound systems 

Three isotopes of the superheavy compound nucleus with Z = 110 were formed by means of 
the following fusion reactions 



58 M + 208p b ^ 266 11Q j|J 
64 M + 208p b _^ 272 11Q p 



40 Ca + 232 Th ^ 272 n0 |g 
40 Ar + 238TJ ^ 278 n0 ^ 
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at different excitation energies between 66 MeV and 186 MeV. The prefission neutron multi- 
plicities were determined experimentally and estimated theoretically within our model. Both 
the experimental and theoretical results are given in Table 1. 

Using Langevin type equations to describe the fusion process [^], [|, we obtained the initial 
spin distribution of compound nucleus. As an example, the differential fusion cross section is 
plotted in Fig. 10 for the reaction 40 Ca+ 232 Th at £'; a b=250.4 MeV as a function of the angular 
momentum of the system. The effective fission cross section is marked by bars in the figure. It 
is seen in Fig. 10 that only the lowest angular momenta contribute to the fusion fission process 
of 272 110 due to the fact that the fission barrier of 272 110 vanishes at L>22h. All isotopes of 
the element Z=110 have a rather small fission barrier which disappears rapidly with increasing 
spin (L). This is illustrated in Fig. 11, where the fission barriers of 266 110 corresponding to a 
few L are plotted as functions of the elongation. 

In our calculations we only took into account those trajectories for which the saddle point 
existed, i.e. we neglected the so-called quasifission events. The number of emitted prescission 
particles is mostly governed by the dynamics from the saddle to the scission point, but of course 
the initial conditions (see Section 3.1) play also an important role. In our model we used two 
different sets of Langevin equations. The first set of equations describes the fusion dynamics 
|5|, |9j while the second one, coupled with the Master equations for particle emission, describes 
the fission process [15]. A better solution would be a model in which the whole dynamics 
from fusion to fission as well as the light particles emission would be described by one set of 
equations. Such an approach could immediately solve the problem of the proper setting of the 
initial conditions for the compound nucleus. 

The experimental and theoretical prefission neutron multiplicities for Z =110 compound 
nuclei are written in Table 1 and additionally they are drawn in Fig. 12 as functions of the 
excitation energy of the system. The dashed lines correspond to the results obtained with the 
wall-and-window friction (j ww ) reduced by 50% while the dotted lines stand for the results 
obtained with the standard value of j ww . It is seen that for the lowest excitation energies 
the agreement of theoretical results with the measurements of neutron multiplicities is good, 
while the discrepancy grows with increasing excitation energy. For the most excited compound 
nuclei our theoretical predictions are too large by 1 to 2 units. This result could mean that 
the effective time which the system needs to pass from the saddle to the scission point is too 
long in our model for a very hot compound nucleus. Any of transport coefficients do not 
depend on the temperature. This could be the reason of the observed discrepancy. In the 
near future we intend to perform a new calculation with temperature dependent transport 
coefficients evaluated within the linear response theory |13, [14]]. The temperature dependence 



of the friction parameter could be important, because as one can see in Fig. 12 the reduction of 
the friction parameter by 50% decrease the neutron multiplicities by 1 to 0.5 units depending 
on the excitation energy of compound nucleus. 

In our model the superheavy nucleus with Z=110 is already cold when it reaches the scission 
point after emission of some neutrons and a-particles. So the shell effects begin to play an 
important role for very deformed shapes. The shell effects could lead to the compact fission, 
i.e. splitting of the Z=110 nucleus into two spherical fragments. In our model the shell effects 
are not present but we simulate partially their influence on the number prefission particles by 
counting only those particles which have been evaporated before reaching the elongation of 
nucleus close to the compact scission point (R12/R0 ~ 1.5). 

Contrary to the case of decay of 126 Ba and 188 Pt compound nuclei a hot superheavy nucleus 
with Z=110 exists and fissions at very low angular momenta (L<20h). This range of L corre- 
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sponds to the linear part of the differential fusion cross-section da/dL (see Fig. 10), so we do 
not expect entrance channel effects on multiplicities of prefission neutrons. 

The Coulomb barrier for emission of the charged particles from Z =110 isotopes is rather 
high so the predicted proton and a multiplicities are much smaller than the neutron multiplic- 
ities as one can see in Table 1. 



5 Summary, conclusions and further developments 

In the first part of the present work we described an extension of our model Jl5]| used so far 



for describing the fission dynamics of an excited, rotating and possibly deformed compound 
system. We considered the decay of compound nuclei produced by the fusion of two heavy 
ions. The initial spin distribution was determined by a simple model calculation || |9| which 
took the possible deformation of the two ions into account |3~6fl . Furthermore, we introduced 
the microscopic classical expression for the emission rates of light particles derived in reference 
and showed the link of this expression with the Weisskopf formulation [p2]| . 



In the second part of the paper we presented a detailed comparison of the calculated neutron 
multiplicities with the ones obtained in different experiments involving systems from different 
regions of the mass table. All five systems discussed in the paper ( 126 Ba, 188 Pt and Z=110 
isotopes) are good examples of different fission and particle evaporation mechanisms. 

The light nucleus 126 Ba fissions at very high angular momenta (L ~ 80h). Its saddle point 
is close to the scission point, so the fission dynamics plays a rather minor role in this case. The 
entrance channel effects influence significantly the decay properties of the excited 126 Ba, via 
different initial spin distributions of compound nucleus. Also huge centrifugal forces and big 
deformations corresponding to the saddle point lead in this case to a true competition between 
the neutron and a-particle evaporation. The experimental prefission neutron multiplicities grow 
with the excitation energy of 126 Ba while it is not always the case in our theoretical estimates. 

In the case of 188 Pt the average angular momentum of the fissioning nucleus is around 60h, 
i.e. 20k less than for 126 Ba. Also the way from the saddle point to the scission point is longer 
so that the dynamics from the saddle to the scission plays here an important role. In addition, 
the evaporation of a few particles does not increase so dramatically the height of the fission 
barrier as in the case for 126 Ba. As a consequence the fission process takes place during a longer 
time and the number of evaporated particles is larger. Similarly as in the case of 126 Ba the 
initial spin distribution (the entrance channel effect) influences the multiplicity of pre-fission 
neutrons which is here too small (by « 30%) in comparison with the experimental data. 

A completely different case is the decay of superheavy compound nuclei with Z=110. Here 
the way from the weakly deformed saddle point to the scission point is very long. The fission 
barrier of a hot compound nucleus with Z=110 is very small and it vanishes already at low 
angular momenta (L<22h). The way in which the Z=110 nuclei has been produced does not 
influence their decay properties since such nuclei exist at small angular momenta only, which 
belong to the linear part of the fusion cross-section da/dL. A majority of compound nuclei 
goes to fission with simultaneous emission of neutrons and a-particles. The number of prefission 
neutrons depends in this case on details of the fission dynamics, e.g. the slope of the potential 
and the magnitude of the friction and diffusion parameters. Our model overestimates slightly 
(~ 15%) the experimental number of prefission neutrons here. The neutron emission process 
cools significantly the Z=110 compound nuclei, so the temperature dependence of the friction 
and diffusion parameters could be "visible" in this case. 
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The present formulation of the model and its application can be improved on several points: 



One may introduce a more detailed description of the formation process of the compound 
nuclei. Such a description will involve explicit trajectory calculations of the incident nuclei 
taking into account the nuclear and the Coulomb potentials for different relative initial 
orientations of the deformed nuclei. Calculations of relative potentials have been worked 



out in this framework j3g, |37[ but, to our knowledge, a full fledged dynamical trajectory 



calculation has yet to be done. 

(ii) The full treatment of fusion fission dynamics with simultaneous emission of the light 
particles will release us from a somewhat arbitrary choice of the initial conditions for the 
compound nucleus. Such a project is at present under investigation. 

(iii) We used a classical description of fluctuations and a form of the fluctuation-dissipation 
theorem which are valid for the case of high temperature. In fact, the temperature of the 
systems which we considered above could be low enough for quantum effects (pairing and 
shell effects) to play a non negligible role. This point has been investigated in the recent 



past [12, [L3L [H within the framework of nuclear transport theory pi. It is our aim to 



come back to this point in the near future. 

(iv) Our fission barrier estimates need also further improvements. The present model bases 
on a modified (see Appendix A) liquid drop formula of Myers and Swiatecki [^] and it 
is known that for the cold light nuclei like 126 Ba it overestimates the fission barriers by 
« 10 MeV |4|. 



For reasons of comparison with experiment we shall generalize our model taking asymmetric 
fission events into account. Many compound nuclei generated from mass asymmetric formation 
process channels appear with a high statistical weight. 

Last but not least it would be nice to extend the test of the model to observables other than 
the light particle emission rates. Out-of plane particle angular distributions and correlations 
between emitted light particles would provide further sensitive tests of the validity of our 
approach. For instance, the ratio of neutrons emitted out of the reaction plane and within the 
reaction plane (defined by the beam and the outgoing fission fragments) would enable us to 
determine the deformation of the emitting compound nucleus. 

Furthermore, we stress the importance of a measurement of the emitted neutrons, protons, 
and a-particles in the same system because the emission rates of these particles are reciprocally 
dependent on each other. 
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Appendix A: Transport parameters used in the Langevin equation 



Fission dynamics correlated with prefission particle emission is generally described in a 
phenomenological framework, by means of quantities such as collective temperature-dependent 
potentials, masses, moments of inertia and friction coefficients. 

The fission process is described here in terms of the unique collective variable q = Ru/Ro- 
The fission path parameterized by q is fixed by means of a minimalization of the free energy of 
a nucleus in the (a , a 2 , a 4 , a 6 ) deformation space at a fixed temperature T. The deformation 
parameters are defined as follows pi 



p\ (z) = E Ml) p 2i (-) , -z < z < z , (Al) 
1=0 z o 

where the function p s (z) is the distance in cylindrical coordinates of any point of the surface 
to the symmetry axis. 

The Helmholtz free energy of the nucleus plays the role of the collective potential and can 
be written, in a semiclassical approximation, as 

F (N, Z, q, L, T) = E(N, Z, q, T = 0) - a(N, Z, q) T 2 + E rot (N, Z, q, L) . (A2) 

The first term on the r.h.s. of eq. (A2) is a liquid-drop (or any macroscopic model) type 
energy expressing the semiclassical energy of the deformed cold nuclear system as a function 
of the mass number A = N + Z, and the isospin asymmetry / = (N — Z)/A. The deformation 
dependence of the free energy is taken into account through the shape functions B s and B cou i. 
The rotational energy is calculated explicitly as E rot = L 2 /2J(q) with a g-dependent rigid-body 
moment of inertia. We have used in eq. (A2) the liquid drop parametrization of Myers and 
Swiatecki ISO . 



The effects of excitation are taken into account in eq. (A2) through the level-density param- 
eter a at a given deformation. As the Helmholtz free energy, which is the variational quantity, 
is given by 

F(T) = E{T) -T S{T) (A3) 



with the temperature T as a Langrange multiplier, one can show p0| , f41f that for a liquid-drop 
type system, entropy and excitation energy read 

S = 2aT (AA) 

and 

E* = E(T) - E(T = 0) = a T 2 (A5) 

from which one obtains the last relation in eq. (A2). 

The level density parameter a can be written in the form 

a(N, Z, q) = a v (l + kj 2 ) A + a s (1 + k s I 2 ) A 2 ^B s (q) + a coul Z 2 A' 1 '* B coul (q) , (A6) 

where the parameters ^=0.0533 MeV^ 1 , ^=0.5261, a s =0.1059 MeV" 1 , A; s =2.7192, and a Co «z=0.000458 Me 



are taken from Ref. |32 |. In this expression, as well as in the macroscopic energy in eq. (A2), 
we have neglected the curvature and compression term proportional to A 1 / 3 since its coefficient 
for the cold system is known experimentally to be small |42| |. When evaluating the Coulomb 



term in (A2) and (A6) the Coulomb exchange contribution is neglected. 
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To describe the fission dynamics one also needs to calculate the mass parameter M (q) which 
enters the kinetic term of the equation of motion describing the fission process. In terms of the 

as 



deformation parameter q the mass is given, in the incompressible fluid approximation |26 

20 

M = Tcpo J \p 2 s {z)A\z) + -B\z)]dz , 

-20 

where 



(A7) 



A(z) 



d 



2(1 



and 



B(z 



p 2 s (z) dq 

Id? 

' 2 3z 



pl{z')dz' 



s A+ dpl 



0q 



(At 



(A9) 



and where po=0.17 fm~ 3 is the matter density. 

The friction coefficient associated with the collective coordinate q is calculated in the frame- 
work of the wall and window model [43| , [|4] . The wall contribution, which is the dominant part 
of the friction parameter 7 associated with the fission mode, is given by 



7T 

7 = g P° v 



2-0 



(f) : 



dz . 



(A10) 



where the average velocity v of the nucleons inside the nucleus is at zero temperature defined 
as 

- = -?- = ?A (3 ^)V3 
c mc 4 mc 

with the Fermi momentum pp = hkp = /i(37r 2 po) 1 ^ 3 - 



(All) 
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Appendix B: Comparison between the microscopic semi-classical formulation of 

emission rates with the Weisskopf formula 



In Ref. ||23|| , a formulation of emission rates has been given which is based on the picture 



of a Fermi gas of nucleons at temperature T being the confined to a deformed rotating square 
well. In this appendix we shall show that in the special case of a spherical well, the model 
of Ref. becomes equivalent to Weisskopf 's emission rate formula [^, ||. We consider the 



example of neutron emission. 

According to Ref. [ 2]J , the total number n of neutrons emitted per time unit is given by 

da f d 3 pf n (f ,p)v ± (f ) w$ (u ± (r )) , (Bl) 



n = 
s 

where da is the infinitesimal element of the surface £ at the surface point r . The momentum 
p and the velocity v of a neutron are related by the equation 

p = mv + mw x r, (-S2) 

where r is the position vector, m is the mass of the neutron, and uj is the angular velocity 
of rotational motion of the nucleus. In the case of a spherical nucleus, a collective rotational 
motion is not possible. Consequently, we may put u = without loss of generality. 

The Wigner function f n (fo,P) then describes a gas of fermions confined to a spherical well 
of depth Vo at a temperature T. It is given by the expression 

fn(f,P) = £;•-. n rP ~ vm ( m ) 



h * l + exp^f^-Vb-^ 



T \2m 

where 9(f) is a step function defined to be 1 for position vectors r of points inside the volume 
fl enclosed by the surface £ 

f 1 forfefiincludingE 
I (J otherwise . 

The parameter \x n is the chemical potential of the neutrons. 

The last factor Wq in formula (Bl) is the semi-classical transmission coefficient for a neutron 
hitting the surface E at the point r*o with a velocity v±(fo) perpendicular to the surface. For a 
spherical surface £ the unit vector normal to the surface at the surface point r*o is given by the 
radial unit vector 

Mro) = ^~, (B5) 

Kq 

where Ro := |r | is the radius of the square well potential. 

In Weisskopf's theory of evaporation, the total number of emitted particles of a given sort 
is represented as an integral over the energy e of the emitted particles |22|, ||. Consequently, 
we would like to rewrite the formula (Bl) as an integral over the variable 

e = ^-V (B6) 
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performing at the same time the remaining integrations in a closed form. For this purpose let 
us first introduce polar coordinates (p, 9 P , tp p ) in momentum space choosing the surface unit 
vector ns(r*o) as the polar axis 

= p sin 6 P cos ip p , 
p v = p sin 6 p sin ip p , 
P( = p cos Op . 

In these coordinates the perpendicular velocity v±(f ) is given by 

p cos P 



m 



(57) 



The transmission coefficient Wq can only be unequal to zero, if the velocity fj_(r ) is positive. 
Thus the polar angle P can be restricted to the interval < P < it/ 2. The total neutron yield 
per time unit is thus given by 



n = 



47T 

7^ 



daldpldOp- ^^^ wt ( P -^ 



m 



m 



T \ 2rr 



-Vb— Mr 



n -1 



+ 1 



(B8) 



The surface integral yields for factor 47t/?q. Instead of integration variables p and P , we 
introduce the energy e and the orbital angular momentum / of the neutron impinging on the 
surface 



I = R p sin Op . 



The Jacobian D of the transformation 



D = 



dp 


de p 


de 


de 


dp 


de p 


dl 


dl 



is obtained from the inverse transformation 

' p= ^2m(e + Vq) 



Op = arcsin (^) 



arcsm 



in the form 



D 



y/2m/(e + V ) 
2^/2mR 2 (e + V )-P 



(B9) 
(BIO) 

(511) 



(512) 



(513) 



In formulating the integration limits for the new variables e and I we must take into account 
the condition 



I' 



< 



£ + V , 



[BU) 



2mR\ ~ 2m 

which implies that the total kinetic energy of the neutron can not be smaller than its rotational 
part. 
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The factor p 3 cos 8 p sin 8 P can be written in terms of the new variables as follows 
p 3 cos 6p sin 9 P = 2m(e + Vq 



I 2 I 



'\| 2mR 2 (e + V ) R 
Making use the Eqs. (B13) to (B15) we obtain for the r.h.s. of Eq. (B8): 

y/2mR 2 (e+V ) 



(515) 



2^4^ J de J dlUjfi^y 

— Vq 



(516) 



2(e+Vo) . / 1 _ P l_ . <(e,Z)- 



v /2mJ?g( e +V; ) )-P ' V 2mflg(e+V ) ' ' [ e (e-M«)/T + ij ' 



or 



n=<&-7de- [ dl.l. ( W ^/ T l) ir (517) 
-Vo o L J 

In Eq. (B17), the transmission coefficient Wq(e,1) is obtained from Wq(v±) by expressing 
the argument v± by the variables e and Z. 



j2mR 2 (e + V ) -I 2 

Introducing the dimensionless angular momentum A instead of the variable / 

l = h\, (519) 

the formula (B17) takes the form 

OO A m ax(£) 

dedX 



B 2 ti 

n= I de I d\——, (520) 



-Vo o 
with 

32^ A \ . n.,cl 



d 2 n 4 X-w cl (e,X) 



dedX h [e^-^)/T + i] 

and 



(521) 



A max (e) := ^praRUe + V ) . (522) 

We wish to compare the result (B21) with the partial width T^(E*,A) for the decay of a 
nucleus of excitation energy E* and vanishing total angular momentum (A = 0) by emission of 
a neutron of energy e a and orbital angular momentum A^ as it is obtained from Weisskopf's 
theory (Refs. p|, ||) 

In (B23), p(E*,A) and p^{E* Rl Ar) are the level densities for the mother and daughter nucleus, 
resp., which depend on the excitation energy E*(E R ) and the total angular momentum A(A^) 
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of the mother (daughter) nucleus. The quantity w n (e a , Xp) is the transmission factor for a 
neutron having the energy e a and the orbital angular momentum Xp. The factor 2 (2Xp + 1), 
represents the product of the degeneracy factors of the emitted neutron and of the residual 
daughter nucleus. We note that, in the derivation of Eq. (B23) and also in our model, effects 
of the spin-orbit coupling on the emission probability are neglected. We now use the level 
densities obtained in the Fermi gas model: 

Here J is the rigid body moment of inertia of the nucleus which is assumed to be the same for 
the mother and daughter nucleus, and "a" is the level density parameter. 
The excitation energies E* and E* R are related by the energy conservation 

E* = E* R + e a - fi n . (526) 

The ratio of the level densities is thus given by 

£f = (exp \2^lE* r - 2 v / o5* 



2 

E* 



(527) 



-2 



<f = (exp [2v^ • (y/l-zg± - l)]) • (l - 
Assuming that the ratio (e a — /i n )/E* is much smaller than 1 

E -^ZJhL « 1 (528) 
E* K ' 

and using the relation between the excitation energy E* and the temperature T as given by 
the Fermi gas model 

E* = aT 2 (529) 
we obtain for the ratio of the level densities 

P(E * R) - exp (- £ -^) . (530) 



p(E*) 1 V T 
Using (B30) we obtain for the neutron width (B23) 

Tf - 2{2 \ + 1} e" 2 -** • w n (e a , Xp) . (531) 

This form of Weisskopf 's general formula agrees indeed with our result (B21) if the transmission 
factor w n is calculated classically (w n = Wq 1 ), if the orbital angular momentum Xp is large 
enough so that 

2Xp + 1 w 2Xp (532) 

and if 

exp ( £ -^M » 1 , (533) 
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so that the Fermi occupation factor in (B21) becomes equivalent to the Boltzmann factor 



1 



T 



(534) 



e t + 1 



e 



We note that the approximation (B34) is already used for obtaining the level density formula 
(B25) from the Fermi-gas model. 

We would like to comment that the well-known formula (B23) is already a slight extension 
of the original result of Weisskopf published in Ref . p2j . There, the transmission factor w n is 
expressed in terms of the cross-section for the absorption of a neutron by the daughter nucleus 
using detailed balance. In this form, Weisskopf 's result is of very general validity. The difficulty 
with it is only that detailed balance relates the transmission factor w n to the cross-section for 
absorption of a neutron by a nucleus with excitation energy E* R and angular momentum A^. 
This absorption cross-section can, of course, not be measured. Therefore, one needs a model 
to calculate the transmission factor. 

The considerations in this appendix demonstrate at the same time that it is meaningful to 
replace the purely classical transmission factor Wq 



by a quantum-mechanical barrier penetration factor as it is done in all our numerical calcula- 
tions. 




(535) 
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Table captions: 

1. Multiplicities of the prefission particles emitted by: 126 Ba 188 Pt and 266 > 272 > 278 HO at dif- 
ferent excitation energies. The theoretical estimates for the compound nuclei with Z=110 
are done for two values of the friction force. In the first rows are the data evaluated with 
the friction reduced by 50% while those in the second row correspond to the standard 
wall-and-window friction. 



Figure captions: 

1. The differential fusion (solid lines) and fission (bars) cross sections for the compound 
nucleus 126 Ba for different entrance channel reactions. 

2. The deformation potential V& ss for 126 Ba as a function of the "fission variable" Ri 2 /Ro from 
the ground state deformation up to the scission point. The different curves correspond 
to the different values of the angular momentum L. The temperature of 126 Ba is fixed at 
T=1.6 MeV. 

3. The same as in Fig. 2, but now the different curves correspond to the different tempera- 
tures T. All plots are done for L=30h. The different curves are shifted vertically in order 
to make the relative changes more visible. 

4. The temperature of the compound nucleus 126 Ba at the different excitation energies as a 
function of angular momentum. 

5. The multiplicities of prefission particles as a function of excitation energy of the compound 
nucleus 126 Ba obtained with the entrance channel 28 Si + 98 Mo. 

6. The same as in Fig. 5 but for the reaction 19 F + 107 Ag. 

7. The same as in Fig. 1 but for the compound nucleus 188 Pt. 

8. The same as in Fig. 2 but for the compound nucleus 188 Pt. 

9. The same as in Fig. 3 but for the compound nucleus 188 Pt. 

10. The same as in Fig. 1 but for the compound nucleus 272 110. 

11. The same as in Fig. 2 but for the compound nucleus 266 110 and temperature T—0. 

12. The experimental (points with error bars) and theoretical multiplicities of prefission neu- 
trons as a function of the excitation energy of the compound nuclei with Z =110 obtained 
in four different reactions. 
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Table 1 



CN 


Reaction 


Elab 


E* 


M exp 


6M™p 




M p 








MeV 


MeV 


- 


- 


- 


- 


- 






204.0 


131.7 


2.52 


0.12 


2.29 


0.03 


0.79 




28 Si + 98 Mo 


187.2 


118.5 


2.01 


0.13 


1.71 


0.00 


0.09 


126 Ba 




165.8 


101.4 


1.32 


0.09 


1.83 


0.00 


0.04 






142.8 


84.1 


- 


- 


0.27 


0.04 


0.88 




iyp + iuy Ag 


147.8 


118.5 


1.85 


0.11 


1.99 


0.00 


0.16 






128.0 


101.5 


1.31 


0.17 


1.80 


0.01 


0.06 




34 S + 154 Sm 


202.6 


100.0 


4.5 


0.7 


3.52 


0.00 


0.10 


188p t 




159.8 


66.5 


2.5 


0.7 


1.10 


0.00 


0.00 




ICQ + 172 Yb 


137.6 


99.7 


5.4 


0.6 


3.79 


0.00 


0.06 






513.9 


185.9 


7.83 


0.46 


8.66 


0.26 


1.04 














9.72 


0.23 


0.95 






486.6 


164.6 


7.35 


0.50 


7.65 


0.20 


0.93 














8.56 


0.21 


0.79 


266 110 


58 Ni + 208p b 


461.1 


144.6 


6.17 


0.48 


6.50 


0.16 


0.81 














7.30 


0.17 


0.73 






436.2 


125.1 


4.74 


0.49 


5.45 


0.10 


0.63 














6.21 


0.11 


0.58 






410.1 


104.7 


4.41 


0.41 


4.35 


0.07 


0.48 














4.94 


0.07 


0.48 






377.0 


78.8 


2.94 


0.36 


2.77 


0.03 


0.30 














3.23 


0.03 


0.33 






472.3 


138.3 


5.98 


0.43 


7.49 


0.08 


0.42 














8.25 


0.08 


0.34 






449.9 


121.1 


5.52 


0.38 


6.40 


0.05 


0.30 














7.12 


0.06 


0.28 




64 Ni + 208p b 


423.0 


100.6 


5.13 


0.33 


5.09 


0.02 


0.20 














5.58 


0.04 


0.20 






403.9 


85.9 


3.54 


0.31 


3.94 


0.01 


0.15 














4.49 


0.01 


0.12 


272 110 




377.0 


65.3 


3.03 


0.32 


2.27 


0.00 


0.10 














2.83 


0.00 


0.08 






351.2 


166.3 


8.40 


0.53 


9.23 


0.13 


0.53 














9.95 


0.15 


0.48 




40 Ca + 232 Th 


298.4 


121.3 


5.81 


0.50 


6.40 


0.05 


0.30 














7.12 


0.06 


0.28 






250.4 


80.3 


3.35 


0.34 


3.56 


0.01 


0.14 














4.04 


0.01 


0.12 






299.6 


127.2 


5.78 


0.51 


7.89 


0.03 


0.14 














8.48 


0.04 


0.15 


278 110 


40 Ar + 23 8U 


280.4 


110.7 


4.96 


0.55 


6.69 


0.00 


0.12 














7.23 


0.03 


0.13 






258.0 


91.5 


4.22 


0.44 


5.08 


0.00 


0.09 














5.59 


0.00 


0.09 
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Figure 3: 
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Figure 6: 
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